Spontaneous VT symmetry breaking and quantum phase transitions in dimerized spin 

chains 
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The occurrence of parity-time reversal (VT) symmetry breaking is discussed in a non-Hermitian 
spin chain. The Hermiticity of the model is broken by the presence of an alternating, imaginary, 
transverse magnetic field. A full real spectrum, which occurs if and only if all the eigenvectors are VT 
symmetric, can appear only in presence of dimerization, i.e. only if the hopping amplitudes between 
nearest-neighbor spins assume alternate values along the chain. In order to make a connection 
between such system and the Hermitian world, we study the critical magnetic properties of the 
model and look for the conditions that would allow to observe the same phase diagram in the 
absence of the imaginary field. Such procedure amounts to renormalizing the spin-spin coupling 
amplitudes. 
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Since the paper of Bender and Bocttcherji it is known 
that non-Hermitian Hamiltonians can display a real spec- 
trum if they are invariant under the joint action of parity 
(V) and time reversal (T) symmetry. The parity opera- 
tor V performs spatial reflection and its action consists 
of changing the sign of both position and momentum, 
whereas the anti-linear time-reversal operator T maps 
p in — p and the imaginary part i in — i. Non-Hermitian 
Hamiltonians that are VT symmetric generate a complex 
extension to the Hermitian quantum mechanics The 
axiom that forces the Hamiltonian H to be Hermitian 
is in fact introduced in order to guarantee the existence 
of a stable ground state and unitary evolutions. How- 
ever, these two basic requirements are still satisfied in 
the presence of 'PT-symmetryi^ 

On the other hand, because of the anti-linear char- 
acter of the parity-time reversal operator, the condition 
[H,VT] — is not sufficient to have a 'PT-symmetric 
system. It is in fact possible to observe spontaneous sym- 
metry breaking in some of the eigenstates of H , associ- 
ated to the appearance of complex conjugate eigenvalues. 
The prototypical model introduced in the literature^ is 
that of system obeying the Hamiltonian p 2 — (ix) . Its 
spectrum turns out to be real only if N > 2. As N 
decreases, the number of complex eigenvalues increases, 
and for N —> 1 + , no real roots are detected. For N < 2, 
the VT symmetry is spontaneously broken. 

The search for non-Hermitian quantum models in dis- 
crete systems has led to study tight-binding particle 
models^ and spin chains^— 

The understanding of such systems on physical bases 
seems to be crucial, given that the meaning of non- 
Hermitian Hamiltonians is rather unclear. Mostafazadch 
showed that a necessary and sufficient condition for the 



reality of all the eigenvalues is the existence of an "equiva- 
lent Hermitian counterpart" ji£ that can be built starting 
from the eigenfunctions of the Hermitian conjugate of the 
initial Hamiltonian. 

It is our aim to study spontaneous breaking of the VT 
symmetry in an exactly solvable model of dimerized spin 
chain. The non-Hermitian term introduced here is a stag- 
gered magnetic field affecting the whole chain. As we will 
show, in order to have a region of parameters with a full 
real spectrum, the nearest-neighbor spin-spin interaction 
needs to be staggered as well. The search for the phys- 
ical bases of such non-Hermitian system will be carried 
out without invoking the "equivalent Hermitian counter- 
part" as described before. Following a more physical ap- 
proach, we will study the critical behavior of the system, 
which can experience a (anti-) ferromagnetic quantum 
phase transition driven by a change in the strength of a 
transverse, real, magnetic field. Then, we will look for 
ordinary Hermitian systems displaying the same critical 
behavior. As we will see, the same phase diagram can be 
obtained considering isotropic XX chains, while in the 
presence of anisotropy the mapping is possible with a 
Hermitian chain with non-linear dependence on the in- 
tensity of the transverse magnetic field. 

The action of time reversal and parity in a discrete 
system is summarized as follows. While the time reversal 
operation T is such that TiT = — i, the effect of the par- 
ity on a system of N spins is such that VafV = <J% +1 _ 1: 
being af the ath Pauli matrix (a = x, y, z) acting on the 
Nth spin. 

We start our discussion by considering a nearest- 
neighbor, Hermitian, dimerized chain of an even number 
N of spin 1/2: 
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with <J% + i = erf . This model has been recently employed 
to study the possibility of observing a completely disen- 
tangled ground state i n i 12 The total Hamiltonian we are 
going to study is given by H = Hq + V, where V is the en- 
ergy due to an imaginary staggered transverse magnetic 
field and can be written as 

N 

v = *f]T(-i)V, (2) 

1=1 

being 77 a real number which measures the deviation of 
H from Hcrmiticity. Without loss of generality, h and 77 
will be assumed to be positive. It is clear that V is VT 
symmetric. Then, one should expect to observe such a 
symmetry in all eigenstates of H . As we will see, if 77 is 
strong enough to make complex some of the eigenvalues 
of H, this symmetry is spontaneously broken. 

The total Hamiltonian H is exactly solvable, as shown 
in Ref.O Neglecting boundary conditions (inessential in 
the thermodynamic limit), the solution can be obtained 
through a sequence of operations. The first step is the in- 
troduction of the Jordan- Wigner transformation, defined 

through of = 1 - 2cjq, = Y[ jKl (l - 2cjqJ q, and 

af = Yij<i fl — ^eje; J cj, which maps spins into spin- 
less fermions.— Then, the chain can be divided into two 
sublattices {cai-i = ai and C21 = b{), and two separate 
Fourier transforms (a/ = (N/2)~ 1 ^ 2 ^ fc exp[— i^ffi] 

and k = (A72)- 1/2 £ fe & fc cxpH^]) can be per- 
formed. Finally, a generalized Bogoliubov transforma- 
tion, mixing a?-, Q>_u, &fej allows one to write 

H = E E A «( fc )fc Q e fe , Q -i), (3) 

where £| , (£k,a,±) are fermionic creation (annihila- 
tion) operators, and 

A±(fc) = \J\ k + fi k ±2y/\ k ij, k -vk, (4) 

with 

A fe = h 2 + 7 2 + 72 - 2 7i72 cos 2k, (5) 

Mfe = Ji+ ,4 + 2J 1 J 2 cos2k-r 1 2 , (6) 

Vk = (J172 + J 2 7i) 2 sin2fc. (7) 

As usual in dimer-type chains, the spectrum consists of 
two branches separated by an energy gap. A + (fc) is usu- 
ally referred as optical branch, while A_(fc) is known as 



acoustic branch. A manifestation of the non Hermitian 
nature of H is represented by the absence of orthogo- 
nality between the quasi-particle modes £k + and £fc . 
A spectral decomposition in terms of orthogonal modes 
would be in fact a signature of Hcrmiticity of H. 

In the case of finite-size chains, we should take into ac- 
count that k can run on integers or half-integers multiples 
of tt/N, according to the parity of the number of excita- 
tions of the eigenstate we are considering^ This char- 
acteristic can be exploited to predict the existence of a 
quantum phase transition in the thermodynamic limitiii 

The excitation energies take a simple expression once 
the anisotropy parameters 71 and 72 are set to zero. Un- 
der these hypotheses, A±(k) reduce to 

A±(k) = h±^. (8) 

Given the structure of Eq. (JSJ, it is easy to find the nec- 
essary and sufficient conditions under which the spec- 
trum of H is real for any k. If J1J2 < 0, the function 
fik reaches its minimum for k — and it is negative if 
rj > \Ji + J 2 \, while when J\ J2 > the minimum corre- 
sponds to k = 7r/2 and negative values appear whenever 
77 > I Ji — J2 1 . Then, we can state that all the eigenvalues 
of H are real numbers if and only if 77 < rj c , where the 
critical value is 

ri c = min{|J 1 + J 2 |,|Ji- J 2 |}. (9) 

Beyond this value, the breaking of the VT symmetry is 
observed independently on the size of the system. From 
Eq. (|9|) wc immediately deduce that, without hopping 
dimerization (Jj = J 2 ), the VT symmetry would result 
unavoidably broken and it would be impossible to obtain 
a full real spectrum. In Fig. [TJ we plot the energies of the 
isotropic chain for three different choices of parameters: 
in one of them the spectrum is real, while in the other two 
cases we observe the appearance of forbidden regions for 
the momentum k. These regions build up around or ir/2 
depending on which of the quantities among | J\ + J%\ and 
\Jy — J2I is lower than 77. Note that, in correspondence 
with the beginning of the VT symmetry breaking, the 
two branches A± touch. 

In the presence of anisotropic exchange (71,72 7^ 0), 
the solution for the critical value of rj is more involved. 
It is possible to state that the condition ^ is still nec- 
essary but not sufficient. The bound for the existence of 
a full real spectrum relies on the positivity of XkfJ-k — 
Imaginary values of A±(fc) could in fact depend on (i) 
Afc/Jfc — Vk < or on (ii) Afe + /ife < 0. But, (ii) implies 
(i), and (i) does not imply (ii). The condition (i) can be 
in fact fulfilled also for positive values of /ik, being Afe 
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FIG. 1: (Color online) Band spectrum (in units of rj) of the 
isotropic spin chain (71 = 72 = 0). Black (dark gray) lines 
correspond to A_(fc), while orange (light gray) lines corre- 
spond to A+(fc). The Hamiltonian parameters are fixed as 
follows: (i) (solid lines) Ji = 2, J2 = 0.4, and h = 1; (ii) 
(dashed lines): Ji = 1.6, J2 = 0.8, and h = 1; (iii) (dotted 
lines) Ji = 1.4, J2 = —0.6, and h = 1. According to Eq. (5), 
the spectrum is fully real only in the case (i), while in the case 
(ii) we have rj > | J\ — J2 | and in (iii) we have rj > | Ji + J2 j . 

and J^fc positive by construction, while, for the condition 
(ii), Hk needs to be strongly negative. In Fig. [21 we plot 
the excitation spectrum of H for different values of the 
parameters. Now, also the value of the field h is crucial 
to determine the symmetry breaking. The disappearing 
of the curves coincides with the presence of imaginary 
roots. 
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FIG. 2: (Color online) Band spectrum (in units of 77) of 
the anisotropic chain. Black (dark gray) lines correspond to 
A-(fc), while orange (light gray) lines correspond to A+(fc). 
We show how the value of h determines the breaking of the 
VT symmetry. We chose Ji = 1.1, J2 = 0.1, 71 = 2.4, and 
72 = —0.8. As for the homogeneous transverse field, we fixed 
h — 0.2 for the solid line and h = 1.5 for the dotted line. The 
VT symmetry breaking can be induced by changing the field. 

In order to gain some information about the meaning of 
non-Hcrmitian Hamiltonians, people started to study the 
properties of the equivalent Hermitian counterpart 1 10 i 15 
However, understanding the connection between this 
mathematical tool and the physics outlined by the model 
is still a demanding point. We then follow a different 
way, trying to analyze the critical properties of H, for 



the region of parameters when its spectrum is real, in 
comparison with those of a Hermitian systems spanning 
the same Hilbert space. 

Beyond the VT symmetry, H is known to be invariant 
under the TL 2 group of the rotations by 7r about the z axis, 
since it commutes with the operator S = ^fLi°'i- Due to 
the discrete character of this symmetry, the eigenstates 
of H are classified depending on the corresponding eigen- 
value of S, which can be equal to 1 or —1. A quantum 
phase transition (i.e., a transition driven by the varia- 
tion of an internal parameter of the Hamiltonian which 
manifests itself in a non-analytical change in some prop- 
erties of the ground state)^ takes place whenever the 
transverse field h reaches a value such that, in the ther- 
modynamic limit, the ground-state energies of the two 
parity sectors become degenerate. The signature of the 
breaking of the symmetry S is given by a non- vanishing 
value of the (anti-) ferromagnetic order parameter. A 
simple way to find the critical region of such models is 
based on the count of the number of excitations of each 
of the two symmetry sectors ^b 17 Two critical fields ap- 
pear that allow for negative excitations energies for some 
special k*. To obtain them, we ask to vanish. Then, 
k* = 0,7r/2. In these cases, 

A_(0) = v^ 2 + (7i - 72) 2 - y/(Ji + J2) 2 - v\10) 

a-(tt/2) = v^ 2 + (71 + 7 2 ) 2 - V(Ji - J2) 2 - ran) 

and the corresponding critical fields are 

h%> - v / (^i+^) 2 -r/ 2 -( 7 i-72) 2 , (12) 
h [ c ] = V(Ji-J2) 2 -V 2 -hi+72) 2 - (13) 

Let us assume, without loss of generality, < Hq\ If 
both of them are real, a symmetry breaking mechanism 
allows for (anti-) ferromagnetism for any field within the 

region h(P < h < hffi ■ In the case of isotropic chain 
(71 = 72 = 0), the model belongs to a different class of 
universality and a symmetry breaking mechanism takes 
place without the existence of an order parameter. 

Our investigation about the relationship between our 
model and the conventional Hermitian quantum mechan- 
ics begins with the isotropic model. Let us assume, for 
the sake of clarity, J\ , J2 > and J\> J2 . According to 
Eq. ©, it is possible to define a Hermitian model, whose 
Hamiltonian will be called Hh, with the whole spectrum 
of H, by introducing two new exchange parameters J[ 
and J' 2 such that 

Mfc = J'i + J2 + 2 J{ J'2 cos 2k. (14) 

This can be achieved through J[ = aJ\ and J' 2 = J-z/a, 
where a is a renormalization parameter which diminished 
the difference with the two hopping amplitudes and then 
the dimer character of the chain. There exist two posi- 
tive values of a (a\ and a 2 ) that allow to satisfy Eq. (|14l) . 
They are related by J[(fli) = J^a?) and J [(0,2) = J^i^i)- 
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Depending on the choice between a\ and 0,2, two Hermi- 
tian models can be built that are connected with each 
other by a translation of one step along the chain. Fur- 
thermore, Eq. (|14[) admits also — a\ and — a 2 as possible 
solutions. These roots would map ferromagnets in anti- 
ferromagnets and vice-versa. 

Performing this kind of analysis, the condition © has 
a simple physical interpretation, since for 77 = i] c we 
find J[ — J2 and Hh represents a perfect translation- 
ally invariant chain. Then, if 77 > ry c , it is not possible 
to identify a suitable H^, given that r\ c corresponds to 
the maximum amount of renormalization we can induce 
in the system. In spite of the renormalization of J\ and 
J2, the critical region of the Hh, bounded by J[ — J 2 and 
J[ + J2, has the same width of that of H. 

By releasing the hypothesis of isotropy, i.e., by consid- 
ering 71, 72 7^ 0, the renormalization procedure described 
before is not enough to find a Hermitian model Hh shar- 
ing with H the same spectrum, due to the dependence on 
J of Vk ■ We are however naturally led to modify also the 
anisotropy parameters 7 to preserve Vk from the renor- 
malization, by defining 7J = 71 /a and 7 2 = 072- Due to 
this transformation, is also modified. Then, a map- 
ping can be realized from H(h) to Hn{h') with h! such 
that h 2 + jf + 7I = h' 2 + 7 12 + Y 2 2 . Thus , we have found a 
Hermitian operator with the same spectrum of H subject 



to a transverse field of non-linear strength. 

To resume, we have introduced a non-Hermitian, VT- 
symmetric, dimerized spin chain, and found the condi- 
tions under which this model admits a fully real en- 
ergy spectrum. In an effort to understand the physical 
meaning of this non Hermitian system, we have analyzed 
its critical properties studying the (anti-) ferromagnetic 
quantum phase transition and looked for Hermitian mod- 
els exhibiting the same phase diagram. Considering an 
isotropic chain, the Hermitian counterpart can be ob- 
tained scaling the hopping parameters. In the absence of 
dimerization, the Hermitian counterpart cannot be found 
independently of the strength of the non- Hermitian term. 
Considering an anisotropic chain, also the anisotropy pa- 
rameters need to be renormalizcd. As a consequence, the 
Hermitian Hamiltonian has a non-linear dependence on 
the strength of the transverse magnetic field. 
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